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Abstract 



A thermodynamical minimum principle valid for photon radiation is shown 
to hold for arbitrary geometries. It is successfully extended to neutrinos, in 
the zero mass and chemical potential case, following a parallel development 
of photon and neutrino statistics. This minimum principle stems more from 
that of Planck than that of classical Onsager-Prigogine irreversible thermo- 
dynamics. Its extension from bosons to fermions suggests that it may have a 
still wider validity. 
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I. INTRODUCTION 



Photon radiation has a distinctive quahty in that it interacts httle enough with matter 
that it is typically far from thermal equilibrium even when matter may be close to equi- 
librium. This is the origin of classical radiative transfer [0,0, which represents a link be- 
tween far-from-equilibrium properties of radiation to the near-equilibrium thermodynamics 
of matter, leading to some curious thermodynamical consequences . These in particular 
concern minimum entropy production in a non-equilibrium steady state (NESS). This is a 
property distinct from the classical theory of irreversible thermodynamics as few of the 
requirements for that theory hold: one does not even have a local thermodynamical bilinear 
form of entropy production HJ^. 

Is this distinctive linkage one that may only be found in the domain of bosons? The 
answer is surely not, as neutrinos, which are fermions, are even more tenuously linked to 
matter than photons. However, a beam of neutrinos, unlike one of photons, is not at all 
the norm. In the case of fermions, in contrast, we must turn to the atypical domains of a 
supernova and the early Universe to have macroscopic beams interacting meaningfully at a 
thermodynamical level with matter. These physical situations involve a boundary in space 
or time between neutrino equilibrium and non-equilibrium. If we allow neutrino production 
without effective absorption, ordinary stars provide another example domain. 

However, the question is to what extent the thermodynamical properties that hold for 
photons also hold for neutrinos. That is, is there a similar minimum principle for entropy 
production in the case of neutrino radiation? This paper addresses this question by gen- 
eralizing "radiation" to include any exactly or nearly massless particles whose number is 
not conserved. We consider only the case of matter in local thermal equilibrium (LTE) and 
where macroscopic thermodynamics is valid, not the more general case when LTE breaks 
down and kinetic theory is necessary 

We proceed here with the assumption that the rest mass and chemical potential of 
neutrinos are zero. Clearly a more comprehensive treatment must include the possibility 
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that neither is zero P]. However there are advantages to proceeding in this manner in the first 
instance. By making this choice we ensure that the neutrinos are as much hke the photons as 
possible. The rest mass is not so much a problem here as is the chemical potential. A non- 
zero chemical potential presents a qualitative departure from the thermodynamic properties 
of photons in that it implies an additional independent thermodynamical variable. However, 
what is learned in this case can be a guide to future work. 

We find that neutrinos have an entropy production minimum principle in the steady state 
similar to that of photons, which also manifests itself as a conservation principle for energy. 
Implicit in weak reactions involving neutrinos is the conservation, not only of electric charge 
Q, but of lepton number L and baryon number B as well. Q, 5, and L conservation in 
weak reactions play a non-trivial role, unlike Q in purely electromagnetic processes. These 
quantities are assumed to be exactly conserved in the microscopic sense. But conservation 
in this paper may also have a secondary thermodynamical meaning: absence of sources or 
sinks of these numbers in the form of macroscopic, thermodynamic reservoirs. 

We proceed in a parallel manner between photons and neutrinos in order to highlight 
differences and similarities and to link with the previous work on photons. 



Phase space for both neutrinos and photons is defined by position r, momentum p, 
implying energy e. In the case of photons it is often customary to introduce frequency v 
as a proxy for energy and wavenumber k instead of momentum. The energy of unpolarised 
particles in the phase space volume d^p d?r is 



where n is the mean occupation number for either photons or neutrinos, h is Planck's 
constant. The entropy of that same volume is 



II. GENERAL DEFINITIONS AND RELATIONS 



2ne 



d^p d^r 



(1) 




(2) 
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where k is Boltzmann's constant. The upper signs correspond to neutrinos and the lower to 
photons. 

A momentum-dependent temperature Tp may be introduced for this small phase volume, 
by forming the derivative of entropy with respect to energy, 

;i = -ln(n-^Tl) , (3) 
J p e 

which takes its physical meaning from the time-independent steady state of noninteracting 
quanta. This makes the cell indistinguishable from one that shares the same temperature 
with all other cells. We extend this naturally to neutrinos given the assumptions of the 
paper: the rest mass rrii, and the chemical potential /ii, are both zero. One finds 

— = --^ dx dVt d'^r , (4) 

h'^ \ he J 

given the new variable x such that e = pc, x — e/kTp, and d^p — p^dpdD,, where dD, is an 
element of sohd angle. 
In thermal equilibrium, 

and we may drop the subscript p on the temperature. The energy in the phase volume may 
then be written 

'^i^^L^dxdnd'r . (6) 
{hey (e^±l) ^ ^ 

From this we note that when an integration over x is carried out that the fourth-power 
law follows for both photons and neutrinos. The only difference between the two is in the 
numerical factor of the integral due to the "±" in the denominator in the integrand of the 
X integration. The entropy of the phase volume is treated similarly. After some simple 
manipulations it becomes 



2k{kT) 
{hcf L(e^±i) 



± x"^ ln(l ± e-^)l dx dQ dh . (7) 



This implies the expected third-power behaviour for entropy in equihbrium, for neutrinos 
as well as photons. 
The four integrals: 



/o (e^±l) 16 U5 

and 

f ±.M,>(l±0.. = i^(l)(g) . (9) 

are easily deduced by simple series expansions. Prom these we find the energy per unit 
volume into solid angle dfl, 

'-^{-)T'da . (10) 

16 ^TTC^ 

where a = 2A;*^7r^/(15c^/i^), the Stefan-Boltzmann constant. Similarly for the entropy, 

(-V-^(-)T'dn . (11) 

The vector flux density of energy into solid angle dfl with direction ih is 

^(-)r'md!2 , (12) 

16 ^TT^ 

and for entropy, 

The integrals (8,9) give the canonical fermion factor of 7/8 relative to bosons. We conclude 
that the flux density per solid angle, known variously as the specific intensity or radiance, 
is for energy. 



16 ^TT 

and for entropy, 

15t1/(7 



nr- , (14) 
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Even out of equilibrium, we may relate the specific intensity le, for a given e, to n from 
equations (1) and (4) and the flux density into a solid angle, 

The specific entropy intensity is 



J. 



2ke 
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The fundamental extensive quantities are the specific entropy flux and specific energy 
flux 7g. Note that expression (3) is recovered by forming dJe/dle- 

Although neutrino number is not conserved, lepton number is, and it is thus physically 
important to define a specific number flux for neutrinos corresponding to I^: 

where 

Ie = eN, . (19) 

In the case of photons, number is not interesting, as photon number is not conserved. In 
the case of fermions, some number conservation law always holds (because of the half-integral 
spins). With zero chemical potential, however, the neutrino number is a purely auxiliary 
quantity and depends on the energy flux. Conversely, we could take neutrino number as 
fundamental and energy as derived; in either case, only one variable is independent. In the 
case of non-zero chemical potential, not treated in this paper, neutrino number and energy 
flux become independent variables, with an energy-dependent chemical potential /i^. 



III. ENTROPY PRODUCTION AND ENTROPY FLOWS 

Entropy is an extensive thermodynamic property and can be localised and integrated to 
determine a global amount. Localization is also possible for entropy production itself. This 
result agrees with the general principle of the locality of physical interactions, so long as a 
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thermodynamical picture is valid. It also means that thermodynamics is not restricted to 
a macroscopic box. If we divide space into distinct regions, boundary surfaces are defined; 
entropy can be moved between regions and the notion of entropy flux across a surface follows. 
The entropy production rate can be expressed in a balance or conservation equation, 

e, = |^ + V-.F, (20) 

where is the entropy production rate per unit volume, s is the volume density of entropy, 
and T is the entropy flux density. 

It is only by convention that this equation is called a conservation equation. It is really 
an accounting equation, with no implication of conservation. In fact, this equation allows a 
general statement of the second law of thermodynamics: > 0. 

A somewhat artificial distinction may be made now in equation (20), between radiative 
and matter processes. In the former case we necessarily must always consider full phase 
space, while in the latter we assume a near-equilibrium distribution in energy, nearly identical 
in all directions (so-called "local thermal equihbrium" or LTE) . In that case we find 

es^eT + e], (21) 

where the superscript m denotes non-radiative components which shall be termed "matter." 
7 denotes radiative components. Writing these entropy sources out explicitly we have. 



3, +^ + V.Y. + V.H, (22) 

where is the volume density of entropy in matter, and is the volume density of 
entropy in radiation. and H denote non-radiative and radiative entropy flux densities, 
respectively. 

As radiation we mean here, of course, photons or neutrinos, while other particles take 
the role of being non-radiative and in LTE. A mixture of near-equilibrium (LTE) matter 
and the far-from-equilibrium radiation is a typical one in the Universe. It is the mixture in 
which you are immersed while reading this page: you are warm, and yet you can read these 
words radiatively. 



By using the equations of state and balance equations for extensive variables, we re- 
express the entropy production rate in the form 

es = Ei^-fc^fc + Va, ■ + + V ■ H , (23) 

k 

where the sum is over contributions from extensive variables with index k. is the conjugate 
intensive variable divided by the temperature, is the creation rate of variable k (for 
example, the rate that internal energy is created from the radiation field, nuclear reactions, 
or viscous dissipation). is the flux density of variable k (for example, the flux density of 
internal energy in the case of diffusion). The prime denotes the value in the rest frame of 
the medium. As massless radiation is without a rest frame, the separation of the entropy 
production into these two parts thus turns out to be not at all artificial. 

If we assume in ( p2D a steady radiation field, and integrate over a finite volume V 
bounded by a surface S, with element dS, containing all of the matter, then the overall 
entropy production rate S is 

because matter fluxes must vanish across S. 

If we ignore matter transport processes, equation (23) becomes 

S = / EK^fc} dV+ fu-dS. (25) 

k 

We may arrive at this result, alternatively, by imagining that the process is steady and that 
the entropy change of matter (the first integral) is reversibly drained off to a heat or another 
type of reservoir. That is, for this steady case, we deal only with a subsystem, thermody- 
namically speaking, and so conservation laws do not hold macro scopically (i.e. integrated) 
in the subsystem alone. Of course this has no bearing on microscopic conservation laws. 

It is worth noting that the terms under the first integral of ( p5| ) are all due to processes 
in matter and are not a part of the entropy production for photons or neutrinos. It is a 
common misconception to interpret the radiation heating rate over the temperature, which 
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is a possible term under the first integral, as the entropy production of radiation. It should 
be clear from this construction that Jy eJdV is all accounted for through the second integral 
of ®. 



IV. MINIMUM ENTROPY PRODUCTION 

Equation (pSf ) provides a structure for computing the entropy production rate due to the 
interaction of matter and radiation for many finite bodies locally in equilibrium. The first 
term represents whatever changes are manifest in the entropy of the body while the second 
term accounts for changes in the radiation field itself due to the interaction. 

If photon radiation is impinging on a body of temperature T in a vacuum, the entropy 
production rate is from (25) just 



V 



V-F 

T 



dV+ / H-rfS 

Js 



(26) 



where the volume V is any containing the body, and F is the volume flux density of energy 
radiation. The latter is the integral over all solid angles and energy of /^n, where n is the 
unit vector defining the direction that (16) flows in. If the temperature is (artificially 
held) uniform over the body, then 



1 



F-dS 



H-dS . 



T Js Js 

If the surface area of the body is A, and it emits as a black body, then 



(27) 



T 



F'-dS 



- aT^A - 



W-dS 



4 , 
+ -aT^A , 
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(28) 



where the remaining integrals represent impinging photon radiation, and the superscript i 
denotes an impinging flow only, which is independent of the state of the body. 
It follows that 



dT 



1 

2^2 



F'-dS 



aT^A 



(29) 



or for a minimum 
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aT^A ^ = . (30) 



s 

That is, the entropy production rate is a minimum in the steady state, implying energy 
conservation but for for an arbitrary geometry and impinging field. 

Equations (26) to (30) could represent the entropy production of a (black) planet irra- 
diated by photons originating elsewhere from a star, or could represent a gas cloud under 
similar photon radiation. They could also represent the entropy production of a target of 
matter irradiated by a laser. One can envisage matter in the laboratory as being connected 
to a heat bath, or circumstances where the heat capacity is very high, to justify the steady, 
isothermal construction of the rate of entropy change in matter (first term in (26)). Obvi- 
ously the isothermal construction may be relaxed \T^\ using the equations of Section III, but 
that is beyond the scope of this paper. 



While these equations represent an extension of previous work [p|,pl| to arbitrary matter 
and radiation geometries, the remarkable thing about them is that they represent a minimum 
entropy production principle that is foreign to the classical minimum entropy production 
principle of Prigogine There are no Onsager reciprocity relations, no "linear" empirical 
fiux-force laws, and no meaningful thermodynamical fluxes and forces. In the latter regard, 
there need only be one temperature defined for the entire problem. 

The question of whether this remarkable property is restricted to boson radiation in the 
from of photons has not been put previously, so we now consider the corresponding problem 
for fermions in the form of neutrinos. As in the case of photons we turn to equation (p5[), 
but at this point the differences between photons and neutrinos emerge, not in the second 
(radiation) term, but in the first (matter) term. That is because of the different manner 
in which neutrinos interact with matter. While photons do not conserve their number and 
trivially conserve electric charge, neutrinos are linked to the conservation of charge, lepton 
number, and baryon number through the structure of weak interactions. Here we consider 
only first-generation fermions and only nucleons for hadrons: 

z/ + n — >• p + e~ , (31) 
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and all related reactions. Thus for neutrinos (UB) becomes 



(32) 



Recall that /i^ = is assumed for neutrinos, hg, hn, and fip represent rates of change of 
number densities for electrons, neutrons and protons respectively, each of which is multiplied 
by its corresponding chemical potential. Assuming chemical equilibrium within the matter. 



This, together with an isothermal and black body (neutrino) assumption, leads to 
1 



- ^(tTM) + J^j; {f^n{nn + rip) + fie{ne - Tip)} dV 



(33) 



7 4 

+ - • -aT^A 

8 3 



(34) 



Conservation of baryon number and charge produce 



-aT^A 



W-dS 



7 4 

- ■ -aT^A . 

8 3 



(35) 



Except for the factors of |, this equation is identical to (28). Thus in minimum entropy 
production, the energy balance steady state 
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- -aTM = 



(36) 



must hold as well. Thus we find that this distinctive minimum entropy production result 
is extended to neutrinos, and thereby beyond the limitation to bosons to at least some 
fermions. 



V. LOCAL AND NONLOCAL REGIMES 

Generally, the interactions of neutrinos and photons with matter are most simply viewed 
as purely local. For statistical physics, however, we also need to count momentum states. 
If we make the matter-radiation separation of section 3 and further assume LTE for mat- 
ter, the matter momentum states can be integrated out, leaving the full phase space only 
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for radiation. This condition allows us to avoid a full kinetic theory calculation [0. At 
this juncture, we have a choice of local versus action-at-a-distance representation for the 
radiation. 

While the radiation exists in its own right, in the event that the overall entropy of the 
radiation field is not changing, we need only be concerned with matter-matter interactions 
mediated by radiation, as the radiation terms can be then be integrated out. Radiation 
then becomes merely a special kind of nonlocal heat transport, and E may be represented 
in a multilocal form. This multilocal form is at least true for the first term in (|25| ) even 
when separate changes do take place in the radiation field, such as in conservative scattering 
processes . In all cases, it is the total S that is minimised. 

In the case of quanta in an extended, continuous medium, a common matter-radiation 
LTE is often valid, with a common matter-radiation temperature T{r). This temperature in 
general is not constant in space. LTE holds to extremely high accuracy for photons inside the 
photospheric surface of a star, for example, but not for neutrinos. The entropy production 
associated with the production and transport (diffusion) of photons is 

^^ = J rfr {(l/2)[4acTV(3fi:^)][V(l/T)]2 + £^/r} , (37) 

where is the photon energy production rate density and k,^ is the opacity (inverse mean 
free path) of the matter against photon diffusion |Tl|. Equation (|37D corresponds to the 



second term in equation (pSj), but written as a volume integral of a divergence, up to but 
not including the photosphere. This is in contrast to photon entropy production discussed 
previously in this paper, in that the photons here are virtually in equilibrium with matter 
and so are diffusive, not radiative. 

At the photospheric surface, a single LTE ceases to hold (see below), and the diffusive 
approximation of equation (^) breaks down. Nonetheless the second term of (^) still 
represents the entropy production in the radiation field, but at the photosphere and outside, 
the photons become radiative. The complete photon entropy production of a star of radius 
R (including its photosphere) is 
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S,,sur = ^aTl^iATiR') (38) 

for the photon radiation released into empty space at an ideahsed sharp surface. (Tgur is 
the photospheric surface temperature.) As the volume of integration is increased picks 
up additional contributions to interactions with more matter, for example, with a planet 



Neutrinos emitted by ordinary stars are quite different from photons: as the interior 
temperatures are not high enough for weak interactions to be in LTE, the neutrinos are not 
emitted in anything like a blackbody distribution, and are not subsequently thermalised. 
Their spectra are instead determined almost exactly by the microscopic reaction spectra 
and emerge essentially unaffected by the neutrinos' subsequent travel through stellar matter 
to empty space. If the emitting star does not absorb neutrinos and the receiving Earth-bound 
detector does not emit neutrinos, the total neutrino entropy production is 

= / dV [ d'p 



emitter 



receiver 



dV I rfV ^ , (39) 



where rip (rip') is the neutrino production (absorption) rate density in real and momentum 
space. In a NESS, rip depends only on ep, not on emission direction p. 

In an equilibrated supernova or the early Universe, on the other hand, the neutrinos are 
emitted and absorbed in LTE. The entropy production below the supernova neutrinosphere 
is a function of a single local temperature: 

^. = 1 {(l/2)[7acTV(6/s:,)][V(l/T)]2 + £,/T} , (40) 

like (p?!), with a neutrino mean free path l/ni, and an extra factor of 7/8 in the diffusion part. 
The total Sj, including the neutrinosphere is analogous to (0). An analogous expression 
can be constructed for photon and neutrino entropy production in the early Universe before 
their respective decouplings, but this would require inclusion of general relativity. 

Even if the neutrinos or photons are emitted and absorbed locally as a gas, the system 
in general is not in equilibrium with a single temperature T or T(r). Multiple temperatures 
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can be defined if eacli system component retains its own LTE, with eacli component spec- 
trum tliermal. For example, a plioton gas witli a frequency- dependent temperature may 
interact with matter of temperature T. Then 

is the local specific energy intensity of photons emitted hy the matter. If > T, then 
Je < 0; if < T, then Ie> 0. Thus is always > 0. Stellar atmospheres provide a related 
example. The radiation has a temperature T^(r), while the various chemical species Xi each 
can have their own Ti{r). Thus 

Again is the local specific radiation intensity emitted by the matter. Contributions such 
as (41) or (42) occur in addition to such gradient terms as (37) or (40). 

In a supernova [|13[, the matter and the neutrinos can have different temperatures T 
and T^, and S^, picks up a contribution analogous to (41). The total S again has other 
contributions, e.g., from gradients of T, T^, and T^, such as (37) and (40). 

Neutrinos can interact among themselves by weak neutral currents and change their 
own phase space distribution without any ordinary matter present. The associated entropy 
production is 



dV I de I dQ I de' I dQ' I,,,{r,e,e',n,Q') x 



(43) 



-T,(r,e,^]) mr,e',n') 

which is local in form, /^^'(r, e, e', Q, Q') is the local doubly-specific radiation intensity of the 
neutrinos "shining" on themselves and is proportional to the neutrino-neutrino weak neutral 
current reaction cross section. S^^ vanishes if thermal equilibrium obtains and there is only 
a single temperature, Tp = Tp/ for all p, p', at each point r. 
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VI. SUMMARY AND CONCLUSION 



A non-zero density of entropy production, E, indicates a local process of an irreversible 
nature. It is distinct from the local density of entropy, and it has no direct implications for 
processes elsewhere. For any given process, one may measure local nearness to equilibrium 
by its magnitude, although relative nearness of two different processes may have no meaning 
by this measure. If a single process is free, one expects that the entropy production rate 
will relax toward zero. If some exterior conditions prevent this dynamical relaxation, the 
process will instead simply relax as far as the exterior arrangements permit. 

This dynamical argument anticipates minimum entropy production, at least for suffi- 
ciently simple systems where thermodynamic quantities have meaning. But it also antici- 
pates, as thermodynamics often does, an analogous static construction of the same principle, 
as the relaxation may be imagined to be constrained to quasi-steady conditions in whatever 
state it may be. 

Thus the principle of minimum entropy production found here is not unexpected. How- 
ever in many respects it is very different from the classical minimum entropy production 
principle of Prigogine . It does not depend on Onsager's principle; there are no "linear" 
empirical flux-force laws, nor is there even a meaningful bilinear construction. 

Under the LTE-NESS assumption, systems of both matter |^ and photons |^,^,|TD|JTT| 
exhibit some type of minimum entropy production. A photon-like principle holds, under 
the same assumptions, for neutrinos, as seen in examples given in sections IV and V. These 
examples can be generalised to many local or continuously varying temperatures. 

Conservation laws play identical roles in all three cases, by constraining the microscopic 
interactions of the quanta. Thus, microscopic energy and momentum are always conserved, 
with all the appropriate macroscopic consequences. Because neutrinos and photons are both 
taken here as massless and not conserved in number, the analogy between these two particles 
can be carried through in most aspects. However, neutrinos are fermions, which always have 
some associated conservation law; in this case, lepton number L. Because the weak inter- 
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actions conserve B — L, baryon number B is also conserved. Electromagnetic interactions 
conserve charge, but since photons themselves do not carry charge, this conservation law is 
trivial in radiative transfer, in contrast to the situation for neutrinos, which do carry L. 

The exact masslessness of neutrinos is not proven experimentally |T^, and a logical 
generalization of our results here is to extend the treatment to massive neutrinos. Although 
we have used an electron chemical potential fie in matter, another generalization left open 
is to include a neutrino chemical potential fi,y. These two extensions will be presented in a 
subsequent publication f^. 
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